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Abstract 

We propose a quantum model for the vacuum filled of virtual parti- 
cle pairs. The main originality of this model is to define a density and 
a life-time of the virtual particles. Compared to the usual QED (p, E) 
framework, we add here the (x, t) space time parameters. We show 
how eo and no originate from the polarization and the magnetization 
of these virtual pairs when the vacuum is stressed by an electrostatic or 
magnetostatic field respectively. We obtain numerical values very close 
to the measured values. The exact equalities constraint the free pa- 
rameters of our vacuum model. Then we show that if we simply model 
the propagation of a photon in vacuum as a succession of transient cap- 
tures with virtual pairs, we can derive a finite velocity of the photon 
with a magnitude close to the measured speed of light c. Again this 
is the occasion to adjust better our vacuum model. Since the transit 
time of a photon is a statistical process we expect it to be fluctuating 
and this translates into a fluctuation of c which, if measured, would 
bring another piece of information on the vacuum. 

When submitted to a stress the vacuum may change and this will 
induce a variation in the electromagnetic constants. We show this to be 
the case around a gravitational mass. It gives a physical interpretation 
of a varying vacuum refractive index equivalent to the curved space- 
time in General Relativity The known measurements of the deflection 
of light by a mass, the Shapiro delay and the gravitational redshift do 
bring constraints on the way inertial masses should depend upon the 
vacuum. 

At last some experimental predictions are proposed. 

1 Introduction 

Vacuum is one of the most intriguing concepts in physics. In Quantum Field 
Theory, the quantum vacuum is not empty and is filled with quantized fields 
with non zero ground energy levels. The uncertainty principle in quantum 
physic leads to fluctuations of their expectation energy values which can be 
associated to virtual particles filling the vacuum. The growing role of the 
vacuum in fundamental phenomenon's is witnessed in the Lamb shift [I], the 



variation of the fine structure constant with energy [2] , or the magnetic mo- 
ment of the electron [3] and of the muon [JJ. All these effects are understood 
as due to the vacuum polarization and are well calculated in the framework 
of QED, or even QCD for the magnetic moment of the muon. The vac- 
uum can also be seen as the triggering actor in the spontaneous decay of 
excited atomic states through a virtual photon stimulating the emission[5]. 
In that particular effect experimentalists were able to change the vacuum, 
producing a huge increase or a decrease [H [7] of the spontaneous emission 
rate. Along the same line of thoughts, Casimir predicted that a pressure 
would be present between electrically neutral conducting surfaces [8]. This 
force has been observed in the last decade by several experiments [9] and is 
interpreted as arising from the modification of the zero-point energies of the 
vacuum due to the presence of material boundaries. 

The vacuum can be stressed by several means and the effects of static 
electric or magnetic fields has been theorized [10] and should reflect, for 
instance, in a modification of the speed of light c. What is done is the cal- 
culation of the departure from c due to these vacuum modifications. The 
value of c is assumed and is not calculated. Not only c but the vacuum 
permittivity eo, and the vacuum permeability /j,q are assumed to be fun- 
damental constants and their experimental values have no known physical 
origin. The goal of this paper is to propose a mechanism arising from the 
quantum vacuum which would lead naturally to these three electromagnetic 
constants. While we were writing up this paper, we found an article recently 
published in jTT] where the authors propose a similar approach to explain 
the physical origin of eo and /Uo- Although the mechanism that they pro- 
pose to derive these two constants is different from the one we will propose 
in this paper, the original idea is the same: the existence of the physical 
electromagnetic constants whose numerical values are simply determined ex- 
perimentally, emerge naturally from the quantum theory. To our knowledge, 
we do not know of any other paper which proposes a direct derivation of 
eo and hq and we do not know any mechanism arising from the quantum 
vacuum which lead to a finite photon velocity c. The most important con- 
sequence in our model is that c, eo and ^o are riot fundamental constants 
but are observable parameters of the quantum vacuum: they can vary if 
the vacuum varies in space or in time. For instance, as shown in section [6l 
the vacuum is modified in the presence of a gravitational field, leading to a 
variation of c, eo and fiQ and a physical interpretation of a varying vacuum 
refractive index equivalent to the curved space-time in a gravitational field. 

The paper is organized as follows. First we describe our model of the 
quantum vacuum filled with virtual charged fermion pairs and we show how 



2 



eo and /^o might originate from the polarization and the magnetization of 
these virtual pairs. Then we show that if we simply model the propagation 
of the photon in vacuum as a series of interactions with virtual pairs, we 
can derive a velocity of the photon with a magnitude surprisingly close 
to the measured speed of light. After that, we study the response of the 
quantum vacuum under the action of a gravitational field and discuss the 
consequences. Finally we present a few experimental tests that could be 
at variance with the standard views and in particular we predict statistical 
fluctuations of the transit time of photons across a vacuum path, which 
translate into fluctuations of c. 

2 An effective description of quantum vacuum 

We propose to model some electromagnetic and gravitational properties of 
quantum vacuum as filled with virtual pairs of elementary particles. We 
consider here only virtual pairs of charged fermions. All species of charged 
fermions are taken into account: the three families of charged leptons e, fi 
and t and the three families of quarks (u, d), (c, s) and (t, b), including their 
three color states. 

We will take sums over the fermion type noted i. We use the notation 
Qi = Qi/ e i where qi is the modulus of the i-kind fermion electric charge and 
e the modulus of the electron charge. 

Virtual pairs of type i are assumed to be produced from virtual photons 
filling the vacuum. Their density is assumed to be limited by the Pauli 
principle. This implies that when a pair decays back into virtual photons, 
the infinite vacuum photon density allows to create immediately a new pair 
at the same place. 

We assume that first order properties can be safely deduced assuming 
all pairs are created with an average energy £j, not taking into account a 
full probability density of the pair kinetic energy. The average energy of the 
pair is taken proportional to its rest mass energy: 



where mi is the rest mass of the fermion i and Kg an unknown constant, 
assumed to be independent from the fermion type, which we keep as a free 
parameter, greater than unity. 

The life-time Tj is given by the Heisenberg uncertainty principle: 



Si = Kg 2m jC 2 



(1) 



~ 2£~ 
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The largest pair life-time is the e + e~ one, with a value in the 1CP 22 s 
range. 

We also assume that the vacuum density iVj of the i virtual pairs is fixed 
by the minimum distance 5i allowed between two virtual fermion in the 
same spin state. This distance is physically limited by the Pauli exclusion 
principle and is of the order of half the Compton wavelength Xf of the 
fermion. We set: 

2 2rriiC 

where the K$ factor is assumed to be the same for all fermion types and is 
expected to be of the order of unity. Due to angular momentum conserva- 
tion, the fermion pairs are coming from the vacuum in the singlet spin state. 
So, there are two fermion-antifermion spin combinations (namely up-down 
or down-up) per cell, and the i-pair density is: 

N i = l (4) 

With Kg = 1, Si ranges from 200 fm for e + e~ pairs to .6 am for tt pairs, 
and Ni ranges from N e ~ 2.8 10 38 pairs/m 3 to Nt ~ 1.5 10 50 pairs/m 3 . 



3 Derivation of the vacuum permittivity 

Consider a parallel-plate capacitor with a gas inside. When the pressure 
of the gas decreases the capacitance decreases too until there is no more 
molecules in between the plates. The strange thing is that the capacitance is 
not zero when we hit the vacuum. In fact the capacitance has a very sizeable 
value as if the vacuum were a usual material body. The dielectric constant of 
a medium is coming from the existence of opposite electric charges that can 
be separated under the influence of an applied electric field E. Furthermore 
the charge separation stays finite because they are bound in a molecule. 
These opposite translations result in no charge in the volume of the dielectric 
and opposite charge appearing on the dielectric surfaces in regard of the 
metallic plates. This leads to a decrease of the effective charge, which implies 
a decrease of the voltage across the dielectric slab and finally to an increase 
of the capacitance. 

In our model of the vacuum the virtual pairs are the pairs of opposite 
charges and the separation stays finite because the electric field acts only 
during the life time of the pairs. In an absolute empty vacuum the induced 
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charges would be null because there would be no charges to be separated 
and the capacitance of our parallel-plate capacitor would go to zero when 
we would remove all molecules of the gas. 

We will see in this section that introducing our vacuum filled by virtual 
fermions will cause its electric charges to be separated and to appear at the 
level of eo ~ 9 pCb/m (~ 5.10 7 electron charges /m) under an electric stress 
E = 1 V/m. 

We assume that every fermion-antifermion virtual pair of the i-kind bears 
a mean electric dipole di given by: 

di = eQiK d 6i (5) 

where 5i the distance between the pairs given by eq. and an arbitrary 
factor assumed universal, which is expected to be of the order of unity. 

If no external electric field is present, 5% points randomly in any direction 
and the resulting average field produced by the collection of these virtual 
dipoles is zero. 

We propose to give a physical interpretation of the observed vacuum 
permittivity eq as originating from the mean polarization of these virtual 
fermions pairs in presence of an external electric field E. This polarization 
would show up due to the dipole life-time dependence on the electrostatic 
coupling energy of the dipole to the field. At distances large compared to 
Si, in a field homogeneous at 5i scale, this energy is di ■ E. From ([2]), one 
gets: 



h/2 h 

~ Si -di-E ~ 2£ i (l-7/ i cos0) ( } 

- — % dE 
with 9 = (di, E) and rji = — — 

Since it costs less energy to produce such an elementary dipole aligned 
with the field, this configuration lasts a bit longer than the others, leading to 
an average dipole different from zero. This average dipole Di is aligned with 
E. Its value is obtained by integration over 9 with a weight proportional to 
the pair life-time: 

Jq di cos 9 Ti(9) 2-7rsin6' d9 
1 = f£n(6) 27rsin# d9 

To first order in E, one gets: 

Di = d^ = ^-E (7) 
3 3£j y ' 
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We estimate the permittivity eo,j due to z-type fermions using the relation 

Pi = e Q ,iE (8) 
where the polarization Pj is equal to the dipole density: 

Pi = NiDi (9) 

Thus: 

= N 1 D 1= 2$_ lKj^Ql 
°> l E 3£i5f 3K 8 m l c 2 5 l { } 

So: 

2K 2 p 2 

z - d n 2 
£CM " 3K £ K 6 he ^ l 

which does depend only on the fermion electric charge, and no more on its 

mass. 

Each species of fermions increases the induced polarization and therefore 
the vacuum permittivity. By summing over all the fermions, one gets the 
estimation e*o of eo: 

K 2 2e 2 



eo 



K £ K 5 3hc 



Ea 2 (ii) 



The sum is taken over all fermion types. We do not have to sum on the 
anti- fermions which are already included in the dipole of the virtual pair nor 
on the two spin states which are accounted for in the density in Equation @ . 
For one generation Q 2 = 1 + 3 x (4/9 + 1/9) = 8/3 (the factor 3 stands 
for color). Each generation contributes. Hence, for the three families of the 
standard model J2i Qi = 8> and one obtains: 

£ °"Wic (12) 

If « Kg ~ Ks ~ 1, one gets eo ~ 4.3 10 -12 F.m -1 , a numerical value 
very close to the observed value eo = 8.85 10~ 12 F.m . If we require an 
exact equality, we obtain a first relation constraining three parameters of 
our model of vacuum 

K% 3 

d - -2.04 (13) 



K^Ks 6Aira 

As in [11], the vacuum permittivity does not depend upon the fermion 
masses, but only upon the number of species and the parameters of the 
proposed model of vacuum. This is at variance with the common idea that 
the energy density of the vacuum is the dominant factor [10] . 
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4 Derivation of the vacuum permeability 



The vacuum acts as a highly paramagnetic substance. When a torus of 
a material is energised through a winding carrying a current I, there is a 
resulting magnetic flux density B which is expressed as: 

B = n nI + [i M (14) 

where n is the number of turns per unit of length, nl is the magnetic intensity 
in A/m, M is the corresponding magnetization induced in the material and 
is the sum of the induced magnetic moments divided by the corresponding 
volume. In an experiment where the current I is kept a constant and where 
we lower the quantity of matter in the torus, B decreases. As we remove all 
matter, B gets to a non zero value: B = fiQnl showing experimentally that 
the vacuum is paramagnetic with a vacuum permeability /j,q = 4ir 10~ 7 N/A 2 . 

We propose to give a physical interpretation to the observed vacuum 
permeability as originating from the magnetization of the charged virtual 
fermions pairs under a magnetic stress, following the same procedure as in 
the former section. 

Each charged virtual fermion carries a magnetic moment proportional 
to the Bohr magneton: 

eQiti 

* = ^7 (15) 

(with the notations defined in section |2J) . 

Since the total spin of the pair is zero, and since fermion and antifermion 
have opposite charges, each pair carries twice the magnetic moment of one 
fermion. The coupling energy of a i-kind pair to an external magnetic field 
B is then — 2/2j • B. 

In a similar way as in the former section (Eq. EJ), the pair life-time is: 

n{0) = h 'l n (16) 
with 8 = (fli, B) 

As in the electrostatic case, pairs with a dipole moment aligned with the field 
last a bit longer than anti- aligned pairs. This leads to a non zero average 
magnetic moment < Mi > for the pair, aligned with the field and given, to 
first order in B, by: 

4u 2 

<Mi>=^B (17) 
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The contribution /2o,i of the i-type fermions to the vacuum permeability is 
given by B = flo^Mi or: 



J_ = Mi 

where Mj is the volume magnetic moment: 

Mi = Ni < Mi > 

This leads to: 

1 8/i? 



(18) 



A*o,t 3£i5f 

Replacing <5j, £j and //j by their expressions, one gets 



(19) 



1 8 ce 2 2 



Qi 



which, again, depends only on the fermion electric charge, and no more on 
its mass. Again this shows that the energy density of the vacuum is not the 
relevant parameter. 

Summing the magnetic moments coming from all fermion types, allows 
to give for fiQ the following estimation: 

ZKjK £ h , . 

A ° = !&flf <20) 
As in the former section, Yli Qi = and one obtains: 

3 h 1 

64 e 2 c 

If Kg ~ Kg ~ 1, one gets £iq ~ 6.4 10~ 7 N.A -2 , corresponding to a numerical 
value very close to the observed value /io = 4-7T 10~ 7 N.A~ 2 . If we require 
an exact equality, we obtain a second relation constraining two parameters 
of our model of vacuum: 

n 2567ra 
K$K £ = —— « 1.96 (22) 



^ = KlK £ -- 2 - (21) 



Using Eq. (j21j) with Eq. (fTTj) from former section, e~o/io is found not to depend 
upon the vacuum energy density but moreover to be independent from the 
number of fermion generations: 

W. = - 2 ^ (23) 
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This implies K d Kg = 2, which means that the pair electric dipole should 
be equal to the fermion charge multiplied by the fermion Compton length, 
whatever the total energy of the pair and whatever the pair density. 



Equations (|13p and (j22[) allow to compute Kg and K d from Kg: 

K.= ™ (24, 



and 

K d = 1.6 ZfK~ £ (25) 

The bound Kg > 1, which requires that the fermion pairs are produced 
with some kinetic energy, converts into bounds on Kg and K d : 

Kg < 1.25 and K d > 1.6 

So far, this model is self consistent and proposes a quantum origin to 
electromagnetic constants. 



5 Derivation of the light velocity in vacuum 

We have shown how the virtual particle model of vacuum explains its elec- 
trostatic and magnetostatic properties. The Maxwell equations impose that 
the speed of light c is thus given by the relation c 2 = l/(eo/io)- However 
no physical process has been proposed to explain the real origin of a finite 
velocity for a photon of null inertial mass. We propose in this section a 
mechanism which leads to a finite speed of light arising from the interac- 
tion of the photon with the virtual pairs. The propagation of the photon 
in vacuum is modeled as a series of interactions with virtual pairs present 
in vacuum. When a real photon propagates inside the vacuum, it interacts 
and is temporarily captured by a virtual pair during a time corresponding 
to the life-time r of the virtual pair. As soon as the virtual pair disappears, 
it releases the photon to its initial energy and momentum state. The photon 
continues to propagate with a bare velocity cq which is assumed to be much 
greater than c. Then it interacts again with a virtual pair and so on. The 
delay on the photon propagation produced by these successive interactions 
corresponds to the finite velocity of light. 

The mean free path of the photon between two successive interactions 
is A = (cA)" 1 , where a is the cross-section for the photon capture by a 
virtual fermion and A is the numerical density of the virtual fermion pair. 
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The mean time for a photon to cross a length A is A/cq+t. It corresponds 
to a photon velocity, noted here v, given by 

A 1 

v = n = i 

- + t i + criVr 

co c 

The bare photon velocity cq is assumed to be much greater than c. This 
bare velocity corresponds to a velocity of light in an empty vacuum with no 
virtual particle. It is therefore natural to assume that cq is infinite. In other 
words, time does not flow between two successive interactions. That means 
that time (as defined as the duration to go from one state or one point to 
another one) would physically arise from the life-time of the virtual particles 
and the uncertainty principle. 

We obtain a general expression of the photon velocity v in function of 
three parameters of the vacuum model: 

1 

v = 

ctNt 



The cross-section a for a real photon to interact and to be trapped by 
a virtual pair of fermions during its whole life-time can be estimated from 
the Thomson cross-section aThomson = jd J (A ) 2 where \ c = h/(m e c) is 
the Compton wavelength of the electron. 

For Thomson interaction, the factor a 2 corresponds to the probability 
a that the photon is temporarily absorbed by the real electron times the 
probability a that the real electron releases the photon. However, in the 
case of the interaction of a photon on a virtual fermion, the second a factor 
must be ignored since the photon is released with a probability equal to 1 
as soon as the virtual pair disappears. Thus we write a as: 

ai = K a ^-aQ 2 (\?) 2 (26) 

with the notation given in section [2j where K a , expected to be of the order of 
unity, is assumed not to depend on the nature of the fermion. For instance, 
the cross-section on one particle of a e + e~ virtual pair is a e ~ 100 barns. 
We notice that the cross-section does not depend on the energy of the initial 
photon. It implies that the vacuum is not dispersive as it is experimentally 
observed. 

The mean free path taking now all fermion types into account is given 

by: 

-=y- 
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So: 

Efe OkN k 

The probability to be trapped on the i-kind fermion for a given interac- 
tion is: 

A 

K = A" 

The average photon trap duration r per interaction, is: 



t = ^2 Pin = K^a-iNiTi 



This corresponds to an average photon velocity v given by the general ex- 
pression: 

" = 7 = E^ (27) 

Using equations ([2]) , Q and ([26]) , one obtains the formula for the average 
velocity of a photon in vacuum: 

^2 c (28) 



K CT 32vra E; Q 



We notice that the photon velocity depends only on the electrical charge 
units Qi of the virtual charged fermions present in vacuum but does not 
depend on their masses, nor on the vacuum energy density. Notice this 
is not relying on the Maxwell relation c = l/^/eo/^o- We mean that the 
Maxwell relation is not necessarily always true. 

We do not sum over the two possible spins of fermions. Indeed a photon 
of helicity 1 (-1 respectively) can interact only with a fermion or an an- 
tifermion of helicity —1/2 (+1/2 respectively) to flip temporarily its spin to 
helicity +1/2 (—1/2 respectively). But we take into account the antifermion, 
and again J2i Qi = 8. 

Thus if we sum the contributions from all the charged fermions, the 
velocity of the photon becomes: 



3K$K £ _ KlK £ 



K a 256na K a 
Using equation ([23]) . one gets: 

1 



0.51 c (29) 



c (30) 
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So: 



K a = 1 (31) 

It is remarkable that this simple model leads to a numerical value of the 
speed of light with the correct order of magnitude. 

The average speed of the photon in our medium being c, the photon 
propagates, on average, along the light cone. As such, the speed of the pho- 
ton is independent of the inertial frame as demanded by relativity, despite 
the infinite photon speed between the trapping sequences. 

6 Derivation of the vacuum refractive index pro- 
duced by a gravitational field 

When submitted to a stress the vacuum may change and this will induce a 
change in the constants. We show this to be the case around a gravitational 
mass and we discuss the coherence of that result with General relativity. 

Einstein was the first who noticed in his own words that "the constancy 
of the velocity of light can be maintained only insofar as one restricts one- 
self to spatio-temporal regions of constant gravitational potential^ [12] . and 
later [13]: "the curvature of rays of light can only take place when the veloc- 
ity of propagation of light varies with position. . . We can only conclude that 
the special theory of relativity cannot claim an unlimited domain of valid- 
ity; its result hold only so long as we are able to disregard the influences of 
gravitational fields on the phenomena (e.g. of light)" . The idea that curved 
spacetime in a gravitational field is equivalent to an optical medium with a 
graded vacuum refractive index in a flat space time was first suggested by 
Eddington [T3]. Rosen |15| has shown how to formulate general relativity 
within the framework of a flat metric. It has been also shown (for instance 
by Pauli [16] or by Landau and Lifshitz [T7] ) that the propagation of light de- 

— 1/2 

scribed in general relativity in a static gravitational field as 5 f g 00 dl = 
is equivalent to the Fermat principle: 5 f n ds = 0. Following this anal- 
ogy, Felice |18j . and later Evans, Nandi and Islam [19] show that in many 
metrics of physical interest, the gravitational field can be represented as an 
optical medium with an effective index of refraction which varies in space. 
They show that the propagation of both photons and massive particles are 
governed by a variational principle which involves the index of refraction 
and which assumes the form of Fermat 's principle for the photon and Mau- 
pertuis' principle for massive particles. For instance, in the Schwarzschild 



12 



exterior metric which applies to the spacetime around an electrically neutral 
static spherical mass M, and far from the Schwarzschild radius, the vacuum 
refractive index is given by 

2GM R s , 

™(0 = i + — r = 1 + T (32) 

where G is the gravitational constant, r is the radial distance to the grav- 
itational centre in a flat spacetime, Cqo is the velocity of light in vacuum 
at infinity (vacuum with no gravitational field) and R s = 2GM/c^ X3 is the 
Schwarzschild radius. This formula is similar to the one already proposed 
by Eddington. A different calculation using the same Fermat analogy has 
been also performed recently giving the same result |20j. 

A direct consequence of this approach is that the speed of light c is not 
constant anymore but varies with the radial distance r as: 

<r) = ^ = Co,(l-^) (33) 
n[r) \ r J 

It is shown in [pj5] and |21j that this new formulation leads naturally 
to the phenomenae predicted by general relativity, namely the deflection of 
light rays by a mass, the Shapiro delay and the gravitational redshift. 

To our knowledge, no physical interpretation has ever been proposed to 
explain the origin of that refractive index. We propose here to show that 
our model of interaction of light with the quantum vacuum gives a possible 
physical origin of a varying refractive index in presence of a gravitational 
field. We study the simplest case of a spherical static gravitational field 
produced by a mass M and we give a possible mechanism leading to the 
right expression of the refractive index as given in Equation (I32p . 

As shown in section [5j the speed of light c(r) at any distance r from the 
gravitational mass can be written as c(r) = 1/ (N(r)T(r)a(r)), where N(r), 
r(r) and a{r) are the density of virtual pairs, their life-time and their cross- 
section for photons capture respectively. Infinitely far from the gravitational 
mass M, the speed of light is Coo = 1/ (-/VqoTooCToo). The vacuum refractive 
index n(r) is thus defined by the relation 

/ \ Coo N{r)a(r)T(r) 
n{r) = dr) = Not (34) 

Before calculating how N(r), r(r) and a(r) should vary when the mass 
M stresses the vacuum, we must first make some hypothesis on which fun- 
damental constants are really constant or can vary since we allow c to vary. 
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We assume that the only fundamental constants in nature which do not vary 
in a gravitational field are the Planck constant h, the electric charge e and 
the gravitational constant G. If Ti and e are constant then a appears to be 
also constant in our model. Indeed Equation (fT2|) shows that eo varies as 
e 2 / (he). The product €q.c is thus constant, and a = e 2 / (A-Ke^Tic) is constant. 
However the mass energy mc 2 and the inertial mass m vary. The way the 
mass energy varies must be consistent with the experimental observations of 
the effect of a gravitational potential difference on the apparent energy of an 
atom. This gravitational redshift effect has been initially observed with 57 Fe 
by Pound and Rebka [22] [23] and confirmed later with high precision with a 
hydrogen-maser [24]. A possible physical interpretation of this effect within 
our formalism is the following. The energy hu of the emitted photon is con- 
stant when it is going closer to the Earth, but its velocity c (and therefore 
its wavelength) varies. The atomic energy level, which is proportional to 
the Rydberg constant R y = a 2 /2 x m e c 2 is smaller when the atom is closer 
to the Earth due to the variation of c and consequently of the mass energy. 
To be consistent with the experimental observations, the mass energy mc 2 
must vary as 

m[r)c (r) = m^c^ = m^c^ I 1 - — J (35) 

In other words, the mass energy of a particle at a distance r from a grav- 
itational mass is equal to its mass energy without gravitation field plus 
the gravitational potential. Comparing this relation with Equation (|33p . 
it appears that the mass energy mc 2 varies as c 1//2 and consequently the 
mass m as c -3 / 2 . This behaviour is in agreement with pionneer models de- 
velopped initially by Wilson |25j and later by Dicke [26j who proposed a 
theory where gravitation could be explained as a result of a varying vacuum 
permittivity near matter. In their formalism, the rest energy mc 2 must vary 
as rriocC 2 ^/ ^/Tq. But, as shown in section [3] Equation (fTTj) . our model leads 
to a vacuum permittivity which varies as c _1 . Thus the relation proposed 
by Wilson is in agreement with a mass energy which varies as c 1 / 2 or an 
inertial mass which varies as c~ 3 / 2 . 

The consequence of a varying mass energy is that the life-time r of the 
virtual pairs is modified by a gravitational field. Without field, equation ([2]) 
can be rewritten as: 

In the presence of a static spherical gravitational field, the average energy 
of the virtual pair depends upon the gravitational potential energy (see 
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Equation (|35p ) and its life-time becomes: 



T(r)= h / 2 - h ' 2 (l-^-X 1 (37) 

T[r) K £ 2m(r)c 2 (r) lf e 2m BO cg (l - °° V 2r J ^ 

Far from the Schwarzschild radius, the Equation (|37p may be approximated 
by: 

r(r) - (l + ^) (38) 

It means that the life-time of the virtual pairs is larger when we get closer 
to a gravitational mass M. 

If gravitation acts in the same direction for fermions and antifermions, 
the virtual pair of charged fermions behaves under a gravitational field as 
a monopole. Therefore the virtual pair should move towards the mass M 
during its life-time r. Suppose a unit volume dS.dr where dS is the unit 
area perpendicular to the distance r from the gravitational mass M. The 
mass of virtual pairs contained in the unit volume is N(r)dSdr2m and the 
gravitational force on the unit volume is thus given by: 

CM 

F(r) = — N(r)dSdr2m{r) (39) 

But the virtual pair composed of fermions cannot move freely in space. The 
presence of other virtual fermions in its neighbouring pushes them back due 
to the Pauli exclusion principle which forbids any overlaping of fermionic 
wave functions with the same quantum states. Such an effect corresponds 
to an effective Pauli pressure which opposes the gravitational force. A pres- 
sure is proportional to an energy density, thus the Pauli pressure should 
be proportional to the energy density N(r)mc 2 of the virtual fermions (or 
antifermions): 

PPauli =r ] 2N(r)m(r)c 2 {r) (40) 

where rj is a proportional constant which must be close to 1. 
At equilibrium, the difference of the Pauli pressure forces on each side of 
the unit volume at any distance r is equal to the gravitational force. It 
corresponds to 

CM 

2rj (N(r + dr) - N{r)) m(r)c 2 {r)dS = ^-N(r)2m(r)dSdr (41) 
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We obtain the differential Equation for the virtual pair density: 

dN 1 GM , 

dr (42) 



N(r) rj r 2 c(r) 2 

Far from the Schwarzschild radius, the integration leads to first order in 
R s /r to 

N(r) = N^exp ( +~) « iVoo (l + (43) 



oo 



It means that the density of virtual pairs is larger when we get closer 
to a gravitational mass M. Consequently the occupancy size 5 of a virtual 
pair is smaller near a mass. From Equation (j3J), 5 varies as 



iVoo / V 37] 2r 

We can now calculate how the cross-section a{r) of photon capture by a 
virtual pair varies along r in the gravitation field. From Equations (J3j> and 
(|26p the cross-section is proportional to u(r) ~ a5 2 . However, as discussed 
before, a appears to be constant in a gravitational field. Therefore a(r) 
varies as 5 2 

The combination of Equations (p4|) . ([38l) . (f43l) and (135]) leads to the 
expression of the vacuum refractive index: 

-fr)-fi + ^fi + ^fi-^ («) 



2r / V rj 2r J \ 3r/ 2r 

With 7/ = 1 /3, one obtains the right expression of the vacuum refractive 
index as derived from general relativity and given in Equation (|32|) . 



Since the light velocity is not constant in this formalism but depends on 
the gravitation potential, we should expect an annual modulation of c due 
to the eccentricity of the Earth orbit. The annual variation of the distance 
R of the Earth to the Sun is AR/R ~ 0.03. Thus the expected annual 
variation of c is 

Ac 2GM & AR ^ io 
c Rc^ R ~ biU 



where Mq is the mass of the sun. 
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Experimental test of local Lorentz invariance and possible variation of 
fundamental constants have been recently performed with high sensitivity 
by comparing the resonance frequency u cso of a cryogenic sapphire oscillator 
(CSO) with a reference frequency vh produced by a hydrogen maser. Results 
of a long-term operation over a period of greater than six years [27] show 
no annual modulation with a sensitivity of the order of {y cso — vh)I v h < 
10 -14 /day. This result is in agreement with our model of varying light 
velocity. Indeed the CSO resonant frequency v cso is proportional to the 
ratio of c over the length of the sapphire crystal L: u cso oc c/L. As c varies, 
the Bohr radius Rb = ^ticqTi 2 / '(m e e 2 ) and consequently L vary as c 1 / 2 . Thus 
v cso varies as c 1 / 2 . The reference maser frequency ujj is refered to an atomic 
energy level. Thus uh varies as the Rydberg constant, which varies as c 1 / 2 
as discussed before. Thus the ratio {y cso — Vb)I v h appears constant in our 
model. 

7 Observational predictions 

7.1 c variation with E and B 

This model allows to make a simple prediction on vacuum properties due 
to the e + e~ pair life-time sensitivity to E and B fields. As shown in sec- 
tions [3] and U] the pair life-times depend on the mean polarization induced 
by external fields. Eq. ([6]) and (I16p can be written as: 




(47) 



since the electric dipoles align on E and the spins on B independently. This 
leads to: 




(48) 



which gives, as a function of model parameters 




1 



(49) 



Setting 




K £ 2m 2 c 3 
eQifr 



(50) 



17 



one gets simply: 



<r i >(E,B) = nM — T — T - T7 ^- 2 (51) 



3 \E C J 3 \E< 

The energy of the electric or magnetic dipoles in fields Ef or Ef/c are equal 
to the fermion pair mass (including kinetic) energy. These fields have the 
smallest values for electrons: respectively 

E c e = 2.6 10 18 ^ V/m 

and 

B% = 8.8 W 9 K £ T 

Propagating this life-time change into the light velocity formula ()27|) . and 
taking only the electron effect into account, leads to an effective refraction 
index in presence of static fields: 

3 \Eg) 3 \Bg t 

These effects are expected to show up whatever the orientation of the fields 
with respect to photon direction. Its numerical value depends on Kg and is: 

,„ 6 W- 27 E 2 (MV/m) + 5 Kr 22 B 2 (T) 
n{E, B) — 1 = -2 

This prediction can be tested in high field laboratories. For instance, a 
forthcoming experiment [28] planning to check magneto-electrical properties 
of quantum vacuum as computed from QED [29] and aiming at sensitivities 
in the 10~ 26 range on indexes, should make use of fields as high as 20 MV/m 
and 15 T, for which we predict a (2 1CT 24 + 10" 20 )/K| effect. 



7.2 Transit time fluctuation 

Quantum gravity theories including stochastic fluctuations of the metric or 
compactified dimensions, predict a fluctuation at of the propagation time of 
photons |30j . However observable effects are expected to be too small to be 
experimentally tested. Hogan[3T] also predicts that the non commutative 
geometry at the Planck scale should produce a Planckian noise correspond- 
ing to a spatially coherent uncertainty in transverse rest frame velocity of 
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photon. Hogan and collaborators are developing an interferometer experi- 
ment to test it. Ellis et al. also predict that vacuum should be dispersive 
meaning that photons of different energies would propagate at different ve- 
locities [32]. The observation of gamma ray bursts (GRB) has been used to 
place limits on such a possible variation of c [33J. 

In our model we also expect fluctuations of the speed of light c. Indeed 
in the mechanism proposed here c is due to the effect of successive interac- 
tions and transient captures of the photon with the virtual particles in the 
vacuum. Thus statistical fluctuations of c are expected, due to the statisti- 
cal fluctuations of the number of interactions N s t op of the photon with the 
virtual fermions and the life-time fluctuation of the virtual pairs. 

The propagation time of a photon which crosses a distance L of vacuum 

is 



where t j & is the duration of the k interaction on an 2-kind fermion. Let 
N s top,i be the mean number of such interactions. The variance of t due to 
the statistical fluctuations of N stoP) i is: 



The life-time should also fluctuate. Assuming exponential probability dis- 
tributions of average values Tj, the variance of the sum of the independent 
U,k is: 



t — ^2 ii,fc 





So 




(53) 



Since N stoP: i = L/Ai, one has: 




(54) 



Using equations ([2]), © and (f2"o]h it comes: 




K a 167ra h ^ Qf 
KjK$ 3 ? ^ m 
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The contribution of each fermion is inversely proportional to its mass. There- 
fore the fluctuations of the propagation time are dominated by virtual e + e~ 
pairs. Summing only on the electrons and positrons, 2~ZjQ 2 / m i = V m e 
since only one fermion spin state couples to a given photon helicity, and the 
anti-fermion is accounted for in the density definition. One obtains: 



0i 



aoVI (55) 



with: 



" K^K's 3 m e c z c y K^Kg 



Replacing ^ and ^ by expressions ([23]) an d ([3T]) allows to give <7o as a 
function of i^: 

/a to" 1 / 2 (57) 



Our model predicts a fluctuation of c and a measurement of this effect 
would deliver precious information on the vacuum. 

The experimental way to test fluctuations is to measure a possible time 
broadening of a light pulse travelling a distance L of vacuum. The figure of 
merit to constrain this model is y/L/At where At is the time width of the 
light pulse. The very fact that these statistical fluctuations should go like 
the square root of the distance implies the exciting idea that experiments 
on Earth do compete with astrophysical constraints as we will see now. 

Constraints from astrophysical observations 

The very bright GRB 090510, detected by the Fermi Gamma-ray Space 
Telescope [33], at MeV and GeV energy scale, present short spikes in the 
8 keV - 5 MeV energy range, with the narrowest widths of the order of 
10ms. Observation of optical after glow, a few days later by ground based 
spectroscopic telescopes give a common redshift of z = 0.9. This corresponds 
to a distance, in the standard cosmology, of about 2 10 26 to. Translated into 
our model this is saying that the fluctuation of c has to be smaller than about 
0.7 'fs to" 1 / 2 . It is important to notice that there is no expected dispersion 
of the bursts in the interstellar medium at this energy scale. 

If we move six orders of magnitude down in distances we arrive to kpc 
and pulsars. Short microbursts contained in main pulses from the Crab 
pulsar have been recently observed at the Arecibo Observatory telescope at 
5 GHz [33J. The frequency-dependent delay caused by dispersive propaga- 
tion through the interstellar plasma is corrected using a coherent dispersion 
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removal technique. The mean time width of these microbursts after dedis- 
persion is about 1 fj,s, much larger than the expected broadening caused 
by interstellar scattering. If this new unknown broadening is correlated to 
the emission properties, this implies fluctuations of c smaller than about 
0.2 f s m~ 1 / 2 . 

In these observations of the Crab pulsar, some very sporadic pulses with 
a duration of less than Ins were observed at 9 GHz [35]. This is 3 orders of 
magnitude smaller than the usual pulses. These nanoshots can occasionally 
be extremely intense, exceeding 2MJy, and have an unresolved duration of 
less than 0.4 ns. The light-travel size cSt ~ 12cm. From this the implied 
brightness temperature is 2 10 41 X. Alternatively we might assume the emit- 
ting structure is moving outward with Lorentz factor 7^ ~ 10 2 — 10 3 . In that 
case, the size estimate increases to 10 3 — 10 5 cm, and the brightness temper- 
ature decreases to 10 35 — 10 37 i^. We recall that the Compton temperature 
is W 12 K and that the Planck temperature is 10 32 i^ so the phenomenon, 
if real, would be way beyond known physics. Another way to look at the 
extraordinary physics behind such events is to consider the energy which 
has to be liberated. It is equivalent to the total annihilation of 100 tons of 
matter in a beach ball volume and in less than a nanosecond. We emphasis 
also two features. Firstly, these nanoshots are contained in a single time bin 
(2 ns at 5 GHz and 0.4 ns at 9 GHz) corresponding to a time width less than 
2/VT2 « 0.6ns at 5 GHz and 0.4/VT2 « 0.1ns at 9 GHz, below the expected 
broadening caused by interstellar scattering. Secondly, their frequency dis- 
tribution appear to be almost monoenergetic. However, if these nanoshots 
are really true, the fluctuations of the vacuum would be a thousand times 
smaller. 

Going down a further 10 17 in distances, which is going from the kpc to 
a few hundred meters, we expect fluctuations in the fs range. Therefore 
an experimental setup using femtosecond lasers sent to a multi-pass cavity 
should be able to detect such a phenomenon and help characterizing the 
vacuum. If nanoshots from the Crab pulsar are true, we would have to 
move from femtosecond to attosecond laser techniques. 

7.3 Annual modulation of eo and fi Q 

We have shown in section [6] that the speed of light c varies when photons 
move towards a mass M (eq. |33|) . Thus Equations (jlip and (j20|) show that 
both eo and uq should increase while we go towards a gravitational mass M, 
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as c : 

~ 10—— (\ 2GM \ 

h Coo V rc^ J 

h 1 / 2GM\ 
Coo V ?~<4 J 

An annual modulation of eo or hq is therefore expected with a relative vari- 
ation due to the excentricity of the Earth orbit given by 

Ae _ A/io _ 2GM@ AR 
e Hq Rclc R 

where M & is the mass of the sun, R is the average distance of the Earth to 
the Sunwith a seasonal variation of the order of AR/R ~ 0.03. Thus the 
expected annual modulation is (Aeo)/eo = (A/io)//^o ~ 6 10 -10 . 

The seasonal variation of eo can be illustrated by the variation of the 
capacity C of two metallic plates of surface S separated by a distance d. 
The capacity is given by C = eoS/d. Since c varies, the Bohr radius Rb = 
4ireoh 2 /(m e e 2 ), and consequently the size of the atoms vary as c 1 / 2 . Thus 
the distance d varies as c 1 / 2 and the surface S as c. Then the capacity C 
varies as c _1//2 corresponding to an annual relative modulation of 3 10~ 10 . 



7.4 mumesic atoms 

This model predicts short distance deviations to standard electrostatics at 
scales much larger than for the QED screening effects in the one photon 
exchange process. 

Section [3] is based on eq. ([6]) valid only at large distance compared to the 
dipole size. The dipole size being found to be equal to the fermion Compton 
length, one expects a decrease of eo up to 12.5% for distances small compared 
to Af » 400/?n, since at such scales the e + e pairs do not contribute to the 
polarization. 

In mumesic Hydrogen, the muon Bohr radius is about 250/m and it has 
been observed an increase in the energy of 5% [36]. This result has been 
interpreted as a need for shrinking the proton size, but it corresponds to 
a 2.5% decrease of eo which is a natural consequence of this model which 
is suited to describe bound states. We predict that this effects is higher in 
mumesic higher Z atoms for which measurements are planned. 
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8 Conclusions 



We proposed a quantum virtual particles model of the vacuum including 
the concept of density and life time of virtual particles. Within this frame- 
work, we have shown how eo and (jlq originate from the polarization and the 
magnetization of these virtual pairs when the vacuum is stressed by an elec- 
trostatic or a magnetostatic field respectively. We proposed that the finite 
speed of light could be due to successive transient captures of photon with 
virtual particle present in the vacuum. Our calculated values for eo, hq and 
c are very close to the measured values. The exact equalities constrain the 
free parameters of our vacuum model. In these proposed model, the vac- 
uum may change when submitted to a stress and this will induce a change 
in the electromagnetic constants. We showed this to be the case around 
a mass. The gravitational stress upon the vacuum induces a change of c 
which is reflected into the well known angular deflection of the photons, the 
Shapiro delay and the gravitational redshift. This physical interpretation of 
a varying vacuum refractive index analogy to General Relativity do bring 
constraints on the way inertial masses should depend upon the vacuum. 
This model predicts several effects: 

Stressing the vacuum with an electric or a magnetic field implies a change 
in the speed of light in the reach of soon forthcoming experiments. 

The propagation of a photon being a statistical process we expect a 
fluctuation of the speed of light. It is shown that this could be within the 
grasp of nowadays experimental techniques and we plan to assemble such 
an experiment. 

Mumesic atoms energy levels probe this model at short distance. 
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